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Time-dependent domain wall solutions with infinitesimal thickness are obtained in the 
theory of a scalar field coupled to gravity with the dilaton, i.e. the Jordan-Brans-Dicke 
gravity. The value of the dilaton is determined in terms of the Brans-Dicke parameter uj. 
In particular, the solutions exist for any ui > and as ui — > oo we obtain new solutions in 
general relativity. They have horizons whose sizes depend on uj. 
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In the cosmo logical models for the early universe based on particle physics and general 
relativity (GR), it is known that these models admit topological defects as their classical 



goes through phase transitions, although the chances of actual existence are not necessar- 
ily high. We however anticipate that GR may not be the ultimate theory to describe the 
very early universe for which we need a quantum gravity. Thus we should look for other 
possibilities how these topological defects can be formed in other gravitational theories. 
Besides, the ultimate theory should not only be able to predict the existence of the exist- 
ing objects but also be able to disprove the nonexistence of the nonexisting objects. In 
string theory context there have been various approaches, but the complete results are still 



In this paper we shall take another possible gravitational theory, namely, the Jordan- 
Brans-Dicke (JBD) theory, which is a gravitational theory with the dilaton in the four- 
dimensional space-time. From the string theory's point of view, the JBD theory with a 
specific Brans-Dicke (BD) parameter is a natural effective gravitational theory before the 
dilaton freezes up. Furthermore, a cosmological model can also be built based on the 
JBD theory 0. Thus it is worth while to investigate the existence of topological defects in 
this context. Cosmic string solutions in this theory were previously studied in ref.Q and 
static thin domain wall solutions have been found in ref. || . In this paper we shall look for 
time-dependent thin domain wall solutions in the JBD theory. 

One subtle point in getting a domain wall solution is that solving the equations asymp- 
totically is not good enough. Unless it satisfies proper matching conditions across the wall, 
it is not a domain wall solution of the system. As noted in ref. || , among many asymptotic 
solutions only one with vanishing BD parameter survives as a true solution in the static 
case. Here we shall adopt the same prescription for the matching condition at the wall 
and look for time-dependent solutions. 

Let us consider the action for the JBD theory 



where u is the BD parameter of the theory. In particular the u = — 1 case corresponds to 
the action of the dilaton gravity from string theory 0. For our purpose we take Sm as the 
action for the real scalar field with a double well potential: 



solutions 0. These topological defects in principle can be formed when the universe 



elusive! !§. 




(1) 
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Sm in particular has a TL^ discrete symmetry $ — > — $. Note that in the static case, this 
discrete symmetry is equivalent to the time reversal symmetry. This is why the domain 
wall structure is related to the CP phases []T(J ]TI). Here we shall take $ to be static, but 
the metric tensor can depend on time. In principle one can generalize to the case in which 
$ depends on time. The result however should be the same in the thin wall case if we 
assume the position of the wall does not change. 
For convenience we redefine the variables as 

then the action now becomes 

S = 1(Lg J d ' x ^ R - + S M [e 2 ^g^\. (4) 

Setting <fi = 0, this action reduces to that of a real scalar field coupled to the Einstein 
gravity, in which case a time-dependent thin domain wall solution is known to exist ||12j| . 

In our case although g^ u is a physical metric of the space-time, eq.([|) is very convenient 
for practical purposes. Although </> is the actual dilaton, we will call <p a dilaton too since 
they are related by a simple field redefinition. We shall also assume that <fi does not depend 
on time. 

By varying eq.(|]) with respect to the new metric g^, we obtain 

R»v = d^dv4> + SnG(T^ - \g^T) , (5) 
where the "energy-momentum" tensor is given by 

t. v = -^=^| = -\ [^tg^djbde* - 2e 2 ?%<S>d v <S> + 9lu *w\{& - v 2 ) 2 } ■ 

9 9 (6) 
Note that is not covariantly conserved due to the dilatonic contribution. The physical 
energy-momentum tensor T™ atter + T$ satisfies the gravitational equation of motion for 
JBD theory pf, 

iV - k^uR = 87re 2 * (T™ atter + TjL) , (7) 



where V M T^ atter = and 



- |V 2 ^ (d^dj- \g^d a ^d a 



(8) 



The second term of eq.(|8|) is proportional to u so that it vanishes if u> = 0, but the first 
term is independent of u. 

The field equations for the dilaton </> is 

□ 2 </> = -Lc^v/^-cg^ = -8nGf3T (9) 

and the matter scalar field satisfies 

d^V^ge^g^d,,)® - 2A v / =^e 4/3 ^$($ 2 - rj 2 ) = 0. (10) 

In general, domain wall solutions are obtained in theories where a discrete symmetry 
is spontaneously broken. Note that the action eq.(^) for the matter field $ has a discrete 
symmetry $ — > — $ so that we can look for domain walls, when this symmetry is spon- 
taneously broken. In the case where domain walls have infinitesimal thickness, we can 
approximate the wanted scalar field to behave as 
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rj if z > 0; 
— rj if z < 0. 



Then we are interested in the the domain wall which separates a space of the $ = rj phase 
from a space of the $ = — rj phase. Such an approximation is in fact reasonable for the 
cases where the Compton wavelength of the test particle is much longer than the thickness 
of the wall. 

For the static cases we used the following ansatz for domain wall solutions ||14j| ||15|| [jSj ■ 



ds 2 = A{\z\){-dt 2 + dz 2 ) + B(\z\)(dx 2 + dy 2 ). 

Note that we have required the reflection symmetry between each side of the wall, which 
is an infinite plane perpendicular to the z-direction at z = 0. Now for time-dependent 
domain wall solutions we shall try 

ds 2 = -A(\z\)dt 2 + C(\z\)dz 2 + B(\z\)U(t)(dx 2 + dy 2 ). (11) 

Again we have the reflection symmetry. The time dependent part of the metric is added 
in g xx = g yy . This ansatz in fact is a generalized version of the time-dependent domain 



wall solutions in GR, in which case C = 1 12 



Now we need a prescription to take care of the matching conditions across the wall. 
Strictly speaking \z\ is not analytic at z = 0. However from physics' point of view this 
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merely is an approximation due to the thin wall assumption. One should treat it as a limit 
case of analytic function. One way to introduce a reasonable analytic property for \z\ is 
to use a step function. Thus we can use the following prescription for \z\ to avoid such a 
difficulty §: 

\z\=z[9(z)-9(-z)\, (12) 

where 9(z) is a step function defined by 6(z) = 1 for z > 0, 9(z) = for z < 0. Then 
d z \z\ = [9(z) — 9(—z)] + 2z5(z). If we promise that d z \z\ shall be multiplied with some 
function of z that does not have a pole at z = 0, we can safely use an identification 
d z \z\ = 9{z) — 9{—z). Similarly, d 2 \z\ = 25{z). The reason we try to be careful about 
such analyticity is to check the consistency of the solutions at the wall, which turns out to 
be important to provide interesting constraints on the solutions. Again we would like to 
emphasize that this should be regarded as an approximation. 
Using this ansatz, eq. (|i~l"l) , we find 



where the prime and the dot denote d z and dt respectively. The corresponding field 
equation R tz = leads to 

A = B. (14) 
With this identification the field equations now become 

+ i £ + - I = ~^GXAe^ (* 2 - , 2 ) 2 , (15a) 
A ' C ' |^ =^ + 87rG(e 2 ^' 2 + ^XCe^i^-rj 2 ) 2 ), (156) 

- \^ + = InGXAUe^ ($ 2 - V 2 ) 2 , (15c) 

8nG(3 (^e 2 ^i$ /2 + 2Ae 4 ^($ 2 - r] 2 ) 2 ^ (15d) 

, , -2AA 3 / 2 C 1 / 2 e 4 ^$($ 2 -r ? 2 ). (15e) 

For thin walls we have $ 2 = i] 2 and $' = away from z = so that we shall first 
solve the above equations away from z = 0, then shall check the consistency at the wall. 
Solving eqs.( |15a,~^ ) for z ^ 0, we obtain 

U(t) = e Kt , (16) 
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and from eq. (|15d| ) we obtain 

Ci/2 

^ = ( 17 ) 

where k, a are constants yet to be determined. 

Using these in eqs. (|15a — d ), we can determine 




C - AA ' 2 (18) 
With such C eq. (pT7|) now can be integrated as 

= , (19) 

where we have taken A(0) = 1, C(0) = 1, and accordingly 0(0). 

Since this eq.(|18|) satisfies eqs.( |15a — dj identically away from z = 0, there is no other 
constraint on the function A so that for any A we can obtain an asymptotic solution. This 
is rather intriguing, although some of them might fail to satisfy the matching conditions 
across the wall. Also some proper asymptotic boundary conditions will distinguish them. 
For practical purposes, it is more convenient to choose C and solve for A to check the 
matching conditions across the wall. 

If we assume (7 = 1, then for a = we reproduce the known solutions in the general 



relativity case|T^|, which can be shown easily from eq(|T8|). For a / there are other 
solutions but the indefinite integral to obtain A cannot be performed analytically. 

From now on we shall concentrate on one particularly interesting case of C = A and 
leave the rest as future exercises. For C = A eq. (|15a| ) is precisely the same as eq. ( |15c|) . 
Now eq.(|l8D can be integrated to yield 

A =pe KZ + qe~ KZ , (20) 

where 

Apqn 1 + \a 2 = 0. (21) 

With the required boundary condition, A(0) = 1, i.e. p+q = 1, we can explicitly determine 
p, q as 



p = U 1 -\l 1 + — )i ?=! 1 + v 1 + ^2i- ( 22 ) 



Thus we have 



C = A — B = cosh/dzl — 



2a 2 

1 H ^-sinhdzl 

6K 2 



(23) 



and from eg . flT9|) , we get 



fin 



3 K_ 

2 a 



1 - 



1 _1_ 2q2 I pK|Z 



1 + 



3 K 
2 a 



1 - 



(24) 



k shall be determined in terms of the energy density so that we are left with a free parameter 
a. Thus these are one parameter family of asymptotic solutions. But as we mentioned 
before, not all of them are true domain wall solutions. 

Now let us check the consistency of the solution at the wall. Using the analytic 
property eq.([T2~D we prescribed, For small k eq.( |15o, c{) reduces to 



2a 



k\/1 + ^S(z) = AtiGX ($ 2 - J] 2 )' 



(25) 



which leads to k > 0. Later we shall find out that this small k assumption is indeed related 
to the weak gravitational field approximation in the sense that the energy density is small. 
Similarly, eq.(|156|) leads to 



2n 2 

OK 



t\2 



(26) 



and eq.( |15dl ) reduces to 



i- 2K v 1+ Si^ )=47rG( ^ 2 



(27) 



The consistency of eqs.(p6[)(|27|) identifies the coefficients of the delta- function so that 
we can determine a in terms of f3 as 



a 



3/3k 



1 - 6/3 2 



(28) 



Thus the solution exists only for (3 2 < 1/6, i.e. u > 0. 

Note that k and r\ have mass dimensions and the Newton's gravitational constant G 
has inverse mass square dimension, while A is a dimensionless coupling constant. Using a 
dimensional analysis for a possible thick wall, if we have GXrj 4 ^> ft-/^ and = 0) = 0, 



the LHS of eqs.(|25|)(p6D are effectively comparable to the RHS by smearing out the delta 
function Thus this is a good approximate solution and the condition in fact corresponds 
to the weak gravitational field limit. 

Finally, k can be determined from the "energy" density as follows: Using eqs. (|25|) -(p7|) 
we can compute the "energy- momentum" tensor eq.(|6|) as 

TV = l!_^ (z)diag(l, -1, -1, 0), k > 0. (29) 

As we promised, the small k implies the small energy density. The corresponding physical 
energy-momentum tensor can be computed from eq.([7]). 

Here we would like to call the readers attention to the fact that we have differential 
equations with the Dirac delta-function. Some may find that this is unreasonable be- 
cause after all the Dirac delta-function is not a function but a distribution. But this is not 
completely unreasonable in field theory when we often need to be careful about the analyt- 
icity. The main intention is not to solve the differential equations in question but to check 
the consistency between equations. In this sense this is a sufficiently good approxima- 
tion. In fact one can be more careful about this situation and can introduce distributional 



energy-momentum tensor in terms of delta-function from the beginning [16]. The result 



however is more or less equivalent because we also have derived the distributional "energy- 



momentum" tensor using our prescription eq.(|12|). We can also further clarify the result 
by introducing infinitesimal thickness of the wall and taking approximation around z = 0, 
although we cannot determine the shape of the solution within this thickness exactly. 

Now the physical metric can be obtained by multiplying the conformal factor (see 
eq.(|D) as ds 2 = e 2/3( ^ds 2 so that we obtain 



where 



d^ 2 = e 2 ^A(\z\) (-dt 2 + dz 2 + e K \dx 2 + dy 2 )) , (30) 



V6(3 + (l - V 1 ~ 6 /3 2 ) e " N 



A(b|) = I cosh/dzl sinhftlzl 

1 y/l - 6f3 2 

The value of the dilaton <p can be obtained in terms of (3, which in turn related to the BD 
parameter u in eq.@, as 

= (3(f, + i In 16ttG. (31) 



In particular, at the wall we have 



0(0) 




1 - y/60 

i + VE/3 



+ 



±lnl67rG. 



Note that as u 



oo, i.e. (3 — > 0, we have new solutions in general relativity 



A(\z\) = e~ K ^. 



(32) 



For a given BD parameter u > we have found time-dependent thin domain wall 
solutions in the JBD theory. For ui = 0, it is known that stable static thin domain wall 
solution exists as proven by the author [|J. This in particular also proves that there is no 
thin domain wall solutions for ui < 0, i.e. in the 4-d analogue of the dilaton gravity. 

As in the GR case, the time-dependent domain wall we have obtained also has a 
horizon at \z\ = - tanh -1 yjl — 6/3 2 . As u increases (i.e. j3 — > 0.), the size of the horizon 
also increases. As uj — > (i.e. /3 2 — > 1/6.), the horizon shrinks and eventually disappears. 
Note that the signature of the metric changes behond the horizon. Thus we should regard 
the region within thee horizon also as where our choice of the coordinate system is valid. 
This however is not unusual. As is known eeeeveen in GR case, the metric of time- 
dependent thin walls become singular as \z\ — > oo. 

Although we have shown only one case explicitly, there are enormous amount of possi- 
ble solutions according to eeeq. (|T3|) . At this moment we do not understand why there are 
so many solutions. We however expect that there should be some classification scheme, or 
some arguments to exclude some of them. Also it will be interesting to search for solutions 
which behave well outside the horizon. 

Much work is needed to understand the structure of all other possible solutions and 
what kind of cosmological implications these dilatonic domain wall solutions have. 
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